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ABSTRACT 

Strong iron fluorescence at 6.4 keV and hard- X-ray emissions from giant molecular 
clouds in the Galactic center region have been interpreted as reflections of a past out- 
burst of the Sgr A* supermassive black hole. Careful treatment of multiple interactions 
of photons in a complicated geometry is essential to modeling the reprocessed emissions 
from the dense clouds. We develop a new calculation framework of X-ray reflection 
from molecular clouds based on Monte Carlo simulations for accurate interpretation of 
high-quality observational data. By utilizing this simulation framework, we present the 
first calculations of morphologies and spectra of the reflected X-ray emission for several 
realistic models of Sgr B2, which is the most massive molecular cloud in our Galaxy. 
The morphology of scattered hard X-rays above 20 keV is significantly different from 
that of iron fluorescence due to their large penetrating power into dense regions of the 
cloud, probing the structure of the cloud. High-resolution spectra provide quantitative 
evaluation of the iron line including its Compton shoulder to constrain the mass and the 
chemical composition of the cloud as well as the luminosity of the illuminating source. 
These predictions can be checked in the near future with future X-ray missions such as 
NuStar (hard X-rays) and ASTRO-H (both iron lines and hard X-rays). 

Subject headings: Galaxy: center - ISM: individual objects (Sgr B2) - radiative transfer 
- scattering - X-rays: ISM 



1. Introduction 

Sgr A*, the supermassive black hole with a mass of 4 x 10^ Mq at the center of our Galaxy 



( Schodel et al.ll2002l : lGhez et al.ll2003l : lGillessen et al.ll2009l ). is currently very faint with a luminosity 
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of lO^^-lO^^ erg (jBaganoff et alJboOlh . which is many orders of magnitude lower than that of 



active galactic nuclei (AGNs). Such quiescent activity of the nucleus seems inconsistent with the 
dense e nvironment in the Gal actic center (GO) from which the black hole could accrete sufficient 
matter (|Melia &: Falckd l200lh. Althou gh Sgr A* displays flares whose luminosity increases by 
a factor of ~100 ( Baganoff et al.ll200ll ). the peak luminosity is still 10® times dimmer than the 
Eddington luminosity of the supermassive black hole. Whether Sgr A* has always been so quiescent 
or previously experienced high activity such as low-luminosity AGNs is an important question in 
addressing the evolution and structure of the Galaxy or the accretion mechanism taking place at 
the center of the Galaxy. 

Indications of past a ctivities of Sgr A* ha ve come from X-ray observations of giant molecular 



clouds in the GC region. ISunyaev et al.l (|l993l ) interpreted X-ray emissions associated with molec- 
ular clouds as reflections of a past outburst of Sgr A * and predicte d a co rrelation of 6.4 keV iron 



Ka fluorescent line emission with molecular clouds. iKovama et al.l (jl996l ) discovered with ASCA 
strong 6.4 keV emissions associated with the Sgr B2 giant molecular cloud, which is located at a 
projected distance of 100 pc from Sgr A*, and indicated that the cloud is possibly irradiated by 
a flare of Sgr A* with a luminosity of 2 x 10^^ erg which occurred 300 years ago. Against 
this X-ray reflection nebulae (XRNe ) interpretation, the iron fluorescence could be generate d by 



collisions with cosmic-ray electrons (jYusef-Zadeh et al.l 120071 ) or protons (jDogiel et al.ll2009l ). In 



order to understand the origin of the fluorescence, time variability of the emission is essential since 



In actual fact, detecti ons of time variability of the iron fluo rescent line from Sgr B2 (IKovama et al 
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20inl) in 



a timescale of a few years have provided strong support for the XRN model. 

In the picture of the XRN model, since X-ray photons from the Sgr A* flare are reprocessed 
via the photoelectric effect or scattered in the dense clouds, the observed X-ray emissions from 
the clouds contain important iri f ormat ion about the flare and the molecular clouds. From this 
viewpoint, ISunyaev &: Churazovl (jl998l ) provided theoretical calculations about iron fluorescence 
emissions from an illuminated cloud. Based on a simplified model, they have realized several 
general effects on time evolution of the spectral and morphological characteristics of the reflected 
X-ray emission. They also pointed out that a low-energy tail structure, or a Compton shoulder, of 
the iron line, which will be investigable by future high-resolution imaging spectroscopy, contains 
information on the time elapsed since the outburst. The strong iron fluorescence is the most 
prominent feature of XRNe and has great potential for detailed measurements of giant molecular 
clouds in the GC region. 

In addition to the iron 6.4 keV lines, hard X-rays scattered via Compton scattering above 12 
keV from the clouds also have good diagnostic power. Above 12 keV, the cross section of photoelec- 
tric absorption is smaller than the Thomson cross section for the chemical composition identical 
with the solar abundances. Since the hard X-rays are less affecte d by absorption than so ft X-rays, 
they reflect the structure of dense regions in the clouds directly. iRevnivtsev et al.l (|2004l ) reported 
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hard X-rays up to ~200 keV associated with the Sgr B2 molecular cloud observed with Inter- 
national Gamma-Ray Astrophysics Laboratory (INTEGRAL). Moreover, continuous monitoring of 
Sgr B2 with INTEGRAL for seven years reveal ed fading of the ha rd X-ray emission whose decay 
rate is consistent with the iron line variability (jTerrier et al.ll2010l ). However, its morphology has 
not yet been obtained due to insufficient angular resolutions and high background levels of current 
hard-X-ray observations. 

Future X-ray missions in the next decade will p rovide new observational information on XRNe 
in the GC region. ASTfiO-ff (jTakahashi et al.ll2010l l. which is the next international X-ray satellite 
scheduled for launch in 2014, will simultaneously realize both high-resolution imaging spectroscopy 
with an excellent energy resolution of 7 eV (full width at half maximum or FWHM) at 6.4 keV and 
hard-X-ray focused imaging up to 80 keV. The unprecedented spectral resolution of the soft X-ray 
spectrometer system (SXS) onboard ASTRO-H will enable us to resolv e the fine structure o f the 
iron line complex including the Compton shoulder. The NuStar mission (jHarrison et al.ll2010l ) will 
also open a window of har d-X-ray imaging, schedu led for launch in 2012. In 2020's International 
X-ray Observatory (IXO) (jBookbinder et al.l l2010l ) will be able to observe the molecular clouds 
with both high spectral and high spatial resolution. 

In order to interpret high-quality observations with such future missions, a high-quality model 
of the X-ray emission from the illuminated clouds is also required. Since giant molecular clouds 
in the GC region are very massive (~ 10^ — 10^ -^o) and the Thomson optical depth reaches 
~1, absorption by heavy elements and scattering by hydrogen molecules and helium atoms in a 
complex geor netry should be t r eated with care. For detailed treatment of a complicated geometry 
of the clouds, iMurakami et al.l (j200d ) developed a numerical model of XRNe to constrain physical 
properties of the clouds by comparing the numerical calculations to observations. This approach, 
however, would be insufficient for high-resolution data due to difficulties with multiple interactions. 
Multiple scattering is essential for modeling the Compton shoulder and scattered hard X-rays in 
such dense clouds. A nother approach is adoption of Monte Carlo simulations to solve the problem 
of radiative transfer (^Leahy fc Creig htonlll99.4 ISunvaev fc Churazovl Il998l : iFromerth et al.lboOlh . 
In general, Monte Carlo simulations treat particle tracking by calculating the propagation and 
interactions of a photon in matter. Although the Monte Carlo simulation is suitable for multiple 
interactions, calculations of illuminated clouds have been limited to a problem of a simple geometry 
with homogeneous matter distribution. Therefore, a Monte Carlo simulation that treats both 
multiple interactions and a complicated geometry is necessary for accurate calculation of the XRN 
emissions. 

In this work, we have constructed a new model describing the X-ray emission from the clouds 
based on the Monte Carlo simulation. Our approach appropriately treats multiple interactions of 
photons in a realistic clou d geometry after a tech i iique of X-ray spectral simulation of high-mass X- 
ray binaries developed by lWatanabe et al.l (l2003l . l2006l ). Based on this model, we have investigated 
prospective diagnostics of giant molecular clouds in the GC region with high-resolution spectroscopy 
of an iron fluorescent line at 6.4 keV and hard-X-ray imaging. The simulation we developed provides 



- 4 - 



the first calculation of X-ray images and spectra from an XRN which has a realistic configuration, 
aimed at practical observations. The future X-ray missions will be able to perform the diagnostic 
observations proposed in this paper with the help of our calculation framework. The present paper 
describes design of the simulation ( ^2.ip . simulated geometry models of molecular clouds ( ^2.2p . 
results of the simulation and discussion on the results and observational prospects for the 

future missions (^. 



2. Methodology 
2.1. Monte Carlo Simulation 



In order to calculate X-ray emissions from illuminated clouds, we constructed a new framework 
based on Monte Carlo simulation. This framework, which is called MONACO (standing for MONte 
Carlo simulation for Astrophysics and COsmology), is a general-purpose code applicable to photon- 
tracking simulations of rnany astrophysical o bjects. We utilize the Geant4 simulation toolkit library 
(|Agostinelli et al.l 120031 : I Allison et al.l 120061 ) for tracking particles and building geometry, while 
physical processes are our original implementations. Geant4, which is widely used for experimental 
particle physics and nuclear science, provides sophisticated handling of complicated geometries such 
as radiation detectors. 

Figured] shows a schematic concept of the Monte Carlo simulation. In the simulation, we build 
a cloud geometry system containing physical properties such as chemical composition and number 
density of the cloud; then in each Monte Carlo trial, a photon with energy Eq and direction 
is generated at time to and position xq. The simulation tracks the particle until it escapes from 
the system, recording information about the last interaction of the photon {Ei,fl-j_,ti,xi), which 
an observer regards as emission. A detailed description of the Monte Carlo simulation and its 
framework is given in Appendix Rl 



to an observer 
(escaping) 



emission 
the last interaction 

{Ei,fti,ti,xi) 




X-ray source 
initial condition 



Fig. 1. — Schematic concept of the Monte Carlo simulation. 



We have developed our original implementation of physical processes adaptable to Monte Carlo 
simulation frameworks based on Geant4. In neutral matter such as molecular clouds, relevant inter- 
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actions of X-rays are photoelectric absorption and scattering by bound electrons in hydrogen and 
helium. Scattering by electrons in heavy elements is negligible because of their small abundances. 
While the simulation tracks a photon by Monte Carlo methods, an interaction occurs according to 
its mean free path; then a status of the photon is changed according to differential cross section of 
the interaction. 

As cross section data ofphotoelectric absorption, we adopt data based on the Evaluated Photon 
Data Library 97 (EPDLOTj^, which is provided along with the Geant4 toolkit as the data of Low- 
Energy Electromagnetic processes version 6.9. A K-shell fluorescent photon is generated with a 
probability of a fluorescence yield, following a photoelectric effect. As relevant atomic properties, 
line energies of K-she l l fluo re scence, K-she ll fluo rescence yields, and K/3-to-Ka ratios are taken 
from [Thompson et al.l (j201Cl ). iKrausd (jl979l ). and lErtugral et al.l (j2007l ). respectively. In addition, 
Ka2-to-Kai ratios are flxed to 0.5. 

Scattering by bound electrons can be divided into three channels according to the flnal state 
of the target electron: Rayleigh scattering (elastic), Raman scattering (excitation), and Compton 
scattering (ionization). The difference in the target between a free electron and a bound electron is 
significant for the fine structure of the iron line complex including the Compton shoulder. The shape 
of the Compton shoulder is smeared by an electron binding effect which is known as the Doppler 
broadening effect due to the finite momentum of the electron bound to an atom or molecule. For 
atomic hydrogen, cr oss sections and differential cross sections of the three channels are provided by 
analytical formulae (jSunyaev Churazovlll996l ). For molecular hydrogen, we adopt an approxima- 
tion that scattering by the molecule is identical to that by atomic hydrogen except that the cros s 
section of Rayleigh scattering is enhanced by a factor of two per electron (jSunvaev et al.lll999l ). 
Since the second most abundant element, helium, does not have an alytic expressions of a tomic 
properties, we use numerical calculations of differential cross sections (jVainshtein et al.lll998l ). Ap- 
pendix [B] provides details of the implementation of the scattering processes. 

For verification of our simulation, we performed a benchmark calculation of an illuminated 
cloud. Results of this simulation can be directly compared with Fig. 9 of Sunyaev & Churazov 
(1998) under identical simulation conditions. In the simulation setup, a point X-ray source was 
located at the center of a spherical uniform cloud that has a radial Thomson depth of 0.2 from the 
center to the surface. Figure [2] shows emerging spectra from the cloud in different time intervals 
when a short fiare occurred at t = 0. The initial spectrum was a power law with a photon index 
of 1.8 in an energy range of 1-18 keV. Our results nicely agree with the former calculation on the 
wide-band continuum and the iron Ka line shape including the Compton shoulder for each time 
interval. In addition, the spectra from our calculation present iron K/3 at 7.1 keV and K-shell 
fluorescent lines of other metals, which were neglected in Sunyaev &; Churazov (1998). 



^littp://www-nds. iaea.org/epdl97/ 



Energy [keV] 



Fig. 2. — Spectra emerging from a spherical uniform cloud with a radial Thomson depth of 0.2 
in which a short flare occurred at t = 0. The spectra are extracted in different time intervals of 
<t < R/c, R/c<t < 2R/c, 2R/c < t < m/c, m/c < t < AR/c, and 4R/c < t < 5R/c (from 
top to bottom), where R is the cloud radius. Each subsequent spectrum is multiplied by 0.05 for 
clarity. This simulation setup is identical to the conditions of Fig. 9 of Sunyaev & Churazov (1998). 



2.2. Model of Sgr B2 



We used the simulation framework described in §2.11 to model X-ray reflection from the Sgr 
B2 giant molecular cloud in the GC region. Sgr B2 is one of the brightest XRNe emitting strong 
6.4 keV iron fluorescence, which also shows variability over a timescale of a few years. Moreover, 
since this cloud has been intensively studied in other wavelengths such as radio and infrared bands, 
its structure and physical properties such as velocity fields and density distribution of molecular 
gas are relatively well measured compared to other molecular clouds. Sgr B2 is, therefore, the best 
target for investigating X-ray reflection from molecular clouds. 

Figure [3] shows the geometrical setup of the simulation. A photon generator or X-ray source 
simulating the Sgr A* supermassive black hole is located at the center. The projected distance 
between Sgr B2 and Sgr A* is fixed to 100 p c. Although a trigonometric parallax of Sgr B2 with 



the Very Long Base Array (IReid et al.ll2009l) and X-ray absorption measurements toward the GC 
region with the Suzaku satellite (jRvu et al.ll2009l ) suggested that Sgr B2 is positioned in front of 
the GC, there still remains uncertainty of the position of the cloud along the line of sight. The line- 
of-sight position is important in studying the structure of the Galaxy and also makes a difference 
in the time of the Sgr A* outburst. We assumed three different line-of-sight positions of the cloud; 
y = +100 pc, pc, and —100 pc, to investigate effects of the relative position of the cloud to the 
center. In the present paper, we assume that the distance to the GC from the Earth is 8 kpc. 
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The Sgr B2 complex is the most massive giant molecular cloud in our Galaxy and has a com 
plicated structure containing several active star forming H II regions. As intensiv e observations in 



radio bands have provided m odels of the cloud's structure (|Lis &: Goldsmithlll99d : iHasegawa et al 



1994 : Ide Vicente et al.lll997l ). it can be presumed that the cloud consists of three components: (1) 
very dense cores with diameters of ~ 0.5 pc containing star-forming regions, (2) a dense envelope 
which seems to correspond with an XRN with a number density of hydrogen molecule of = 10^- 
10^ cm~'^, and (3) a large, diffuse component with an almost constant density of rin^ ~ 10'^ cm^'^ 
extending to a diameter of 45 pc. 

For simplicity we constructed a spherically symmetric geometry model composed of a single 
dense core and an envelope that has a diameter of 10 pc and a realistic density distribution. We 
determined the diameter of the envelope according to Hasegawa e t al. (1994) and the v alue well 

). In the 



agrees with X-ray morphology obtained by Chandra observations (IMurakami et al 



2001 



model the radial profile of the number density in the cloud is represented as 



10*^ cm" 



(core: r < 0.25 pc) 



'^o I — ) cm~^ (envelope: 0.25 pc < r < 5.0 pc) 



(1) 



The mass of the core is 2 x lO'' Mq, which is a very small fraction of the total mass of the cloud. 
Table [J shows parameters of our cloud models in the simulation. We prepared four models with 
different no values corresponding to total cloud masses from 2.5 x 10^ Mq to 2x IO^Mq (Models 1-4). 
In M odels 1-4, the power-law index a = 1 was assumed, as rad io observations have reported values 
of 2 \Us fc Goldsmith! ll99ol ^ or 0.87 (Ide Vicente et al.lll997l ). To investigate effects of different 
density profiles, we built a model with a = (Model 5) and a model with a = 2 (Model 6), fixing 
the total mass to 5 X IO^Mq. While we assume d a metal abundance of 1.5 protosolar value in 



Models 1-6 as a standard value in the GC region (jNobukawa et al.ll20ld ). we additionally checked 
two models of different values of 1.0 protosolar (Model 7) an d 2.0 pro t osolar (Model 8). The values 
of the protosolar (solar system) abundances are taken from iLoddersI (|2003l ). We ignored the third 
component surrounding the dense envelope to reduce computation costs; instead we introduced 
absorption of the initial spectrum by the third component as 



F{E) (X exp(-AfHCTabs(^))-E"^ (1 keV <E< 400 keV), 



(2) 



where A'^h = 6 x 10^^ cm~^ is an equivalent hydrogen column density of the surrounding diffuse 
component, and aahs{E) is the photoelectric absorption cross section at energy E per hydrogen. 
The photon index 7 of the initial spectrum is 1.8, which is consistent with flares of Sgr A* or X-ray 
emission from Seyfert galaxies. 
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Fig. 3. — Geometrical setup of the simulation. We assumed three different positions of the Sgr B2 
cloud, fixing the projected distance from Sgr A* to 100 pc. 



Table 1: Parameters of the cloud models: power-law index a of the radial density profile of the 
cloud envelope, the normalization of the power law, the metal abundance of the cloud, and the 
total mass of the cloud calculated by these parameters. 



# 


a 


no 


metal abundance 


-^cloud 






[10"^ cm-3] 


[protosolar] 


[IO^Mq] 


1 


1.0 


2.5 


1.5 


2.5 


2 


1.0 


5.0 


1.5 


5.0 


3 


1.0 


10.0 


1.5 


10.0 


4 


1.0 


20.0 


1.5 


20.0 


5 


0.0 


1.5 


1.5 


5.0 


6 


2.0 


13.1 


1.5 


5.0 


7 


1.0 


5.0 


1.0 


5.0 


8 


1.0 


5.0 


2.0 


5.0 



- 9 - 



3. Results 



3.1. Morphology 

We investigated morphologies of iron 6.4 keV line emissions and hard X-ray emissions by 
extracting images from the simulation results. Since X-ray emissions extend over the whole cloud, 
a duration of the illuminating flare of Sgr A* must be longer than the light crossing time of the 
cloud. As the exact value of the duration does not significantly affect the emissions unless it 
is shorter than the light crossing time, we assumed that the flare began at t = —65 yr, kept 
a constant luminosity, and then ended at t = 0. The duration of 65 years is twice as long as 
the light crossing time of the 10 pc large cloud. The X-ray refiection from the cloud becomes 
brightest just before the rear end of the direct light front reaches the cloud edge confronting the 
illuminating source. Time of the brightest moment depends on the relative position of the cloud to 
the flare source and the observer, as summarized in Table [2l The assumed luminosity of the flare 
is Lx = 1.3 X 10'^^ (d/lOO pc)^ erg in an energy range of 1-10 keV, where d = yx^ + y^ is the 
distance between the cloud and the illuminating source. This is normalized so that the iron line 
flux of Model 2 (y = pc) at the brightest moment agrees with an observational value of 1.0 x lO"'^ 



photons s ^cm ^ obtained with Chandra in 2000 (jinui et al.ll2009l ). 



Table 2: Time of the brightest moment elapsed since the end of the flare (t = 0) for different cloud 
positions. 



position 


time of the brightest moment 


[pc] 


[yr] 


(-100, -100) 


122 


(-100, 0) 


303 


(-100, +100) 


756 



Figure m depicts iron 6.4 keV line morphologies at the brightest moment for the eight different 
cloud models and three different line-of-sight positions. The X-ray illuminates the cloud from the 
right side of each image. The case of the middle position (y = pc) is shown in the middle row of 
Fig. m In this case, the cloud becomes brightest at t = 303 yr. As the cloud mass increases (from 
Model 1 to 4), the cloud becomes opaque, making the morphology semicircular. Corresponding 
hard-X-ray images extracted in an energy range between 20 keV and 60 keV are shown in the middle 
row of Fig. m The hard-X-ray morphologies are obviously different from the iron fluorescence due 
to their large penetrating power. Hard X-rays are therefore useful for probing the structure of the 
dense molecular clouds. Images of Models 5, 2, and 6 show the difference due to the density proflle 
of the cloud. Since the cloud with a steeper radial profile (Model 6) has denser inner layers and 
thinner outer layers, its morphology becomes more concentrated and concave. In Model 5, the 
small opacity of the uniform envelope exposes the central core in the hard X-ray band. The metal 
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abundance in the cloud also changes the iron line emission while the hard- X-ray images are similar. 

If the line-of-sight position of the cloud is different, the incident direction of the illuminating 
photons changes and significantly different morphologies appear. The top rows in Figures S] and [5] 
show X-ray morphologies at the brightest moment, t = 756 yr, in the front-illuminated case or the 
position of y = -1-100 pc. The delay of the brightest moment is due to light propagation via the 
cloud behind the GC. In this case, the iron line emissions form into a round shape since a large 
proportion of the cloud front is illuminated by the flare. On the other hand, the bottom row of 
Figures [Hand [5] show morphologies in the back-illuminated case or the position of y = —100 pc, and 
the reflection becomes the brightest at t = 122 yr. Morphologies of the iron line form into a concave 
shape; moreover, the shape becomes a crescent if the cloud is more massive than ~ 1 x 10^ Mq. In 
contrast to the iron emission, morphologies of scattered hard X-rays do not largely depend on the 
relative position of the cloud or the direction of the illuminating source unless the cloud is Compton 
thick. 



Model 1 



Model 2 



Models 



Model 4 



+100 pc 



pc 



-100 pc 




Model 5 



Model 6 




■6-4-2 2 4 6 



>l > 



Model 7 Model 8 




-2 2 4 



Fig. 4. — Morphologies of iron Ka 6.4 keV line emissions for Models 1-8 at the brightest moment. 
Each row corresponds to a different line-of-sight position of the cloud marked at the left. The colors 
are mapped in linear scale and are common in all the images. The size of each square image is 
12 pcxl2 pc, corresponding to an angular size of 5'. 2 x 5'. 2. 



Since Sgr B2 displays fading of the X-ray emission (IKoyama et al.l l2008l : llnui et al.l l2009l : 



Terrier et al.l l2010l ) , we would observe evolution of the X-ray morphologies in the XRN model 



as shown in Fig. [6l At 6.4 keV the central region of the cloud will become darker in the fading 
phase than outer regions if the cloud is massive enough to absorb most of the incident soft-X-ray 
photons in the central regions. In contrast, hard X-rays do not suffer from such significant absorp- 
tion; thus morphological disagreement between the iron line and the hard X-ray is notable in the 
decay phase. An apparent light front intersecting the cloud from the observe r moves fast when 
the cloud is positioned in front of the flare source (jSunvaev &: Churazovlll998l ). For the cloud at 
y = —100 pc morphological evolution is very fast even in the hard X-ray band (Fig. [6] bottom). 



- 11 - 



Model 1 



Model 2 



Model 3 



Model 4 



+1 00 pc 



pc 



-100 pc 




Model 5 



Model 6 




-6 -4 -3 2 4 6 ^-4 -2 2 4 6 

no 



o 




Model 7 



Models 




Fig. 5. — Same as Fig. [4] but for hard X-ray emissions in an energy range of 20-60 keV. 
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Fig. 6. — Time evolution of the morphology of the iron line (left panels) and the hard X-ray (right 
panels) for Model 2 from the brightest moment at intervals of five years. The observation time is 
marked at the top of each image and the cloud position along the line of sight is marked at the left 
of each row. The colors are mapped in linear scale and are common in all the images of the same 
energy band. 
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3.2. Spectral Properties 

In addition to an image of the reflected emission, its spectrmn contains important information 
about the cloud and the flare. Figure [7] shows spectra extracted from the simulation results for the 
models with different masses (Models 1-4). While the assumed incident spectrum is a power law 
with a photon index of 1.8 absorbed by the outer diffuse region (A'^h = 6 x 10^^ cm~^), absorption 
due to the interstellar medium is not considered here. Since scattered hard X-rays above 20 keV 
do not suffer from absorption significantly, their flux increases with the cloud mass. In contrast 
to the straightforward behavior of hard X-rays, dependence of the iron line on the cloud mass is 
complicated due to heavy absorption in dense clouds. In addition to the iron lines, fluorescence 
lines of other heavy elements (e.g. Cr, Mn, Co, Ni) can be observed with high energy resolution. 
Spectral evolutions are considerably different in each cloud mass since the absorption plays an 
important role in the fading phase of the reflected emission. A more massive cloud displays faster 
decay, enhancing the line equivalent width and the Compton shoulder. 

In order to evaluate the spectra quantitatively, we extracted the hard-X-ray flux (20-60 keV), 
the iron line flux, the equivalent width of the iron line, and the shoulder-to-peak ratio of the 
Compton shoulder as a function of time. The iron line was divided into a peak and a shoulder. 
The peak flux and shoulder flux are integrated over energy ranges between 6.0 keV and 6.38 keV, 
and in a range between 6.38 keV and 6.42 keV, respectively. A continuum is evaluated by fitting 
to a power law over a range between 5.5 keV and 6.8 keV without several fluorescent line energies, 
and is then subtracted from both fluxes. The equivalent width and the shoulder-to-peak ratio were 
calculated by using these values obtained. 

Figure [8] shows the spectral parameters as a function of time elapsed since the end of the 
Sgr A* flare. After t = 303 yr, the hard-X-ray and iron line fluxes gradually decrease while the 
equivalent width and shoulder-to-peak ratio increase with time. Although the iron line flux is 
almost independent of the cloud mass in the brightest phase, its decay speed depends on the mass 
because of absorption. The equivalent width of the iron line keeps constant and does not depend 
on the cloud mass when the whole cloud is illuminated. In the fading phase, however, its variation 
significantly depends on the mass; it shows a rapid incre ase in the dense clouds. On e reason for 



this is a difference in photoelectric absorption probability (jSunyaev Sz Churazovlll998l ). In a dense 
cloud, continuum photons around 6 keV are largely and rapidly absorbed though hard X-rays 
above 10 keV, which are responsible for producing the iron line, are not so much. As a result, 
a large equivalent width of the iron line appears. Since the Compton shoulder is generated via 
scattering of the iron fluorescence, it is also a good probe of the Thomson depth of the cloud or the 
total mass of the cloud. The equivalent width and the shoulder-to-peak ratio also depend on the 
metal abundance in the cloud, as shown in Fig. [9j The equivalent width is nearly proportional to 
the metal abundance. A large metal abundance increases photoabsorption probability, suppressing 
scatterings to generate the shoulder. Note that the metal abundance controls amounts of all metals, 
not only iron; changing only the iron abundance could bring different results. 
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Fig. 7. — Spectra of the reflected emissions for the cloud masses of 2.5 X lO^M© (Model 1; black), 
5.0 X IO^Mq (Model 2; red), 1.0 x lO^M© (Model 3; blue), and 2.0 x lO^M© (Model 4; green). 
The thick solid lines are spectra at t = 303 yr (the brightest moment) and the thin solid lines are 
spectra sd, t = 318 yr (fading phase). The cloud position is at y = pc. Left: -wide band spectra. 
Right: enlarged around the iron complex -with a bin -width of 5 eV. Several fluorescent lines are 
seen: Mn Kai(5.90 keV), Cr K/3 (5.95 keV), Fe Kai (6.40 keV), Mn K/3 (6.49 keV), Co Kai (6.93 
keV), Fe K/3 (7.06 keV), and Ni Kai (7.48 keV). 
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Fig. 8. — Hard X-ray flux in a range of 20-60 keV (top), iron line flux (second), equivalent width in 
unit of keV (third), and shoulder-to-peak ratio (bottom) of the Fe Ka line as a function of time for 
cloud masses of 2.5 x IO^Mq (Model 1; black), 5.0 x IO^Mq (Model 2; red), 1.0 x lO^M© (Model 
3; blue), and 2.0 x W^Mq (Model 4; green). The flux unit is 10~^ photons s~^cm~^. The cloud 
position is at y = pc. 
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Fig. 9. — Same as Fig. [8] but for different metal abundances of 1.0 (Model 7; solid), 1.5 (Model 2; 
dashed), and 2.0 (Model 8; dotted). 
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For different line-of-siglit positions of the cloud (Fig. [TOl llip . the dependences of the iron hne 
flux on the cloud mass are different from the cloud at y = pc, while other parameters behave 
in the same way. Since the flux dependence is a monotonic function of the cloud mass, these 
behaviors of the iron line are somewhat easier to interpret than those of the cloud at y = pc. 
For the front-illuminated cloud at y = +100 pc, it is an increasing function because it depends 
simply on the total mass of matter that generates fluorescence. However, the line flux is saturated 
at Mcioud 1 X 10^ Mq. For the back-illuminated cloud at y = —100 pc, on the other hand, 
the dependence is a decreasing function due to heavy absorption in a dense cloud. Moreover, the 
line-of-sight position makes a difference to the decay time of the hard- X-ray emission without any 
significant biases of the cloud structure. In contrast, the time evolution of the iron line flux depends 
not only on the position but also strongly on the amount of absorbing matter. The decay time 
of the hard X-rays to decrease by half from the brightest moment is 13, 23, or 28 years for the 
different positions along the line of sight of y = —100, 0, or +100 pc, respectively. 




Fig. 10. — Same as Fig. [8] but for the cloud position of y = +100 pc. 
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Fig. 11. — Same as Fig. [8] but for the cloud position oi y = —100 pc. 
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4. Discussion 

In the previous section, we describe images and spectra of X-ray reflection from Sgr B2 obtained 
by the simulation. By comparing observations with the simulation results, we can in principle 
constrain physical properties of both the illuminated clouds and the illuminating flare of Sgr A*. 
Since there is no available data of high-resolution spectra and hard- X-ray morphologies, it is beyond 
the scope of this paper to obtain observational constraints. Future data with a high signal-to-noise 
ratio would require more simulation runs with tuned model parameters. Furthermore, it would be 
necessary to assume more realistic geometries that would not have spherical symmetry. In these 
cases, our simulation framework allows us to build a more likely geometry of the cloud and to 
calculate emissions from it. 

Although the XRN model predicts that X-ray emissions from Sgr B2 will get very dim in 10-20 
years, it is fortunate that we are able to observe the response of the cloud to the incident X-rays 
from Sgr A* as a time variability of its morphology and spectrum. The response provides us with 
information about the structure and position of the cloud, as described in the previous section 



The concave shape of the iron line emission from Sgr B2 obtained with Chandra (jMurakami et al 



2001 



) suggests that the cloud position along the line of sight is not behind the GC. Moreover, the 



fast decay time of ~8 yr obtained with the INTEGRAL hard X-ray obse rvations ([Terrier et al 



20101 ). which seems comparable with that of the iron line (jlnui et al.ll2009l ). also supports a near- 
side position of the cloud. This is consistent with the radio trigonometric parallax and the X-ray 
absorption measurements (mentioned in ^2.2p . A more precise estimate will be obtained by future 
measurements of the apparent velocity of the rear end of the light front, which is determined 
by the relative position of the cloud to the observer and the illuminating source. The degree of 
absorption estimated by the slight distortion of the observed morphology of the iron line suggests 
that the mass of the cloud's envelop e is smaller than the value of ~ 1 x W^Mq estimated by radio 



lope 

3). 



observations (iHasegawa et al.l 119941 ). This estimation could be inaccurate if the cloud had a more 



complex geometry such as a non-spherical or highly clumpy structures. 

It is important to check the feasibility of the observations discussed in this paper with fu- 
ture instruments. Here we describe brief prospects of the hard- X-ray imaging and high-resolution 
spectroscopy for ASTRO-H. The diameter of the Sgr B2 envelope which we investigated with sim- 
ulations, 10 pc, corresponds to 4'. 3 in an apparent angular size. An angular resolution of 1.7' 
(half-power diameter or HPD) of the hard X-ray imaging system (HXI) onboard ASTRO-H allows 
us to pick up a few points in a spatial profile of hard-X-ray brightness across the cloud. Despite 
the inadequate resolution to trace the fine structure of the cloud, the resolution is sufficient to 
constrain the total cloud mass and position. Moreover, NuStar will have an angular resolution of 
43" (HPD) or 7". 5 (FWHM) to resolve the cloud structure in the hard- X-ray band. 

For high-resolution spectroscopy, we present observation simulations of Sgr B2 with the SXS 
onboard ASTRO-H foi an exposure time of 200 ksec in Fig. [T2j To demonstrate the advantage of 
the SXS over conventional observations, the same simulations with the Suzaku observatory, which 
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has the current best spectral performance for XRNe studies, are superposed. These predictions are 
based on the cloud of Model 2 positioned at y = pc at the brightest moment {t = 303 yr) and 
15 years later (t = 318 yr). Since the angular size of the Sgr B2 cloud may be larger than the 
field of view of the SXS (3' x 3'), we assumed that the cloud has a brightness averaged over the 
cloud size of 10 pc and has uniform brightness distribution in the entire field of view for the SXS 
simulations. The energy resolution of the SXS is sufficient to resolve the fine structure of the iron 
Ka line including its Compton shoulder at 6.4 keV, iron K/3 at 7.1 keV, and fluorescent lines of 
less abundant metals, overcoming the soft detector response of current X-ray instruments. HM The 
exposure time of 200 ksec is adequate to evaluate the line properties and continuum levels, though 
it depends on the actual flux of the Sgr B2 emission which will decrease. In real observations of 
Sgr B2, the spectrum should be contamin ated by hot plasma emissions of a temperature of 6.5 



keV (jKovama et al.ll2008l : iRvu et al.ll2009l ) composed of a continuum, Hea lines of FeXXV at 6.7 
keV, and Lya lines of FeXXVI at 6.9 keV. For accurate estimation of the contamination level, 
the excellent energy resolution of the SXS will be essential to measuring the plasma temperature, 
which determines the continuum emission via bremsstrahlung. This estimate of the hot plasma 
contamination will be applicable to X-ray observations with current instruments to obtain more 
precise values of the equivalent width of the iron line when we safely assume the invariable flux of 
the thermal emission. 

In addition to the Sgr B2 complex, other molecular clouds in the GC regions are also valuable 
XRNe for studying the Sgr A* activity and the structure of the molecular clouds. The fact that a 
part of the molecular clouds around Sgr A display time variability of the iron fl uorescence provide s 



strong evidence of the past activity of Sgr A* and information on its history (jPonti et al.ll2010l ). 
Our model of XRNe enables us to constrain the positions and physical properties such as total 
masses of clouds in the GC regions. Based on the constraints thus obtained, we could construct 
a consistent view of the structure of the GC molecular zone and the history of the activity of the 
Sgr A* supermassive black hole from a few thousand years ago to the recent past several decades. 
Moreover, there are giant molecular clouds that are faint in X-ray near Sgr A*. Future missions 
with a large effective area will be able to investigate such clouds by means of a Compton shoulder 
which is a remnant of past illumination. 



5. Conclusions 

We have developed a new calculation framework of X-ray reflection from giant molecular 
clouds in the GC region (e.g. Sgr B2) based on Monte Carlo simulations. Our model carefully 
treats radiative transfer including multiple interactions in the complicated geometry of a cloud. By 
utilizing this simulation framework, we investigated X-ray diagnostics of XRNe aimed at hard- X-ray 
imaging and high-resolution spectroscopy for which future missions such as ASTRO-H will open a 
window. Scattered hard X-rays from the cloud show signiflcantly different morphologies from the 
iron fluorescence due to their large penetrating power into a dense cloud. Thus, hard-X-ray imaging 
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Fig. 12. — Observation simulations of Sgr B2 based on Model 2 (y = pc) with ASTRO-H SXS 
(red points) and Suzaku XIS 0+3 (black points) for an exposure time of 200 ks. The left panel 
shows spectra at t = 303 yr (the brightest moment) and the right panel at i = 318 yr (fading 
phase). The simulated data are represented in points with error bars showing la statistical errors, 
superposed on solid lines representing the model convolved with the detector response. The data 
are re-binned for clarity. Besides iron Ka at 6.4 keV, iron K/3 7.1 keV and nickel Ka at 7.5 keV 
are clearly detected. 

will provide new important information on the structure of the illuminating molecular clouds. High- 
resolution spectroscopy also brings a quantitative approach to constraining the physical properties 
of the cloud such as the total mass and the chemical compositions of metals, as well as the luminosity 
of the illuminating source. A large-scale view of X-ray reflection from molecular clouds in the GC 
region will reveal the structure of the GC molecular zone and the history of the activity of the 
central black hole in our Galaxy. 

We thank Engine Churazov for providing us with data concerning differential cross sections 
of scattering by He atom. H. Odaka is supported by research fellowships of the Japan Society for 
the Promotion of Science for Young Scientists. This work was supported in part by Global COE 
Program (Global Center of Excellence for Physical Sciences Frontier), MEXT, Japan. 



If we assumed that radiation did not affect the surrounding medium. X-ray reprocessing (flu- 
orescence or scattering) in the system would be regarded as a linear response of the system. The 
linear nature of the system greatly simplifies the problem because the final results are obtained 
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by convolving the response function, which describes how the system changes the initial photon 
into the final state of the photon. In other words, the problem is reduced to obtaining a Green's 
function of a linear operator which describes the system. X-ray reflection by cold matter can be 
treated in this manner. 

Thus, the purpose of the Monte Carlo simulation is to calculate a conditional probability 
density of emitting a photon with an energy Ei and direction fli at time ti and position Xi under 
initial conditions of a photon emitting with an energy Eq and direction CIq at time to and position 
xq. Figure [T] shows a schematic concept of the simulation. In each Monte Carlo trial, we track 
a photon from certain initial conditions until the photon escapes from the system, recording the 
initial condition and the conditions of photons at the last interaction before escaping. After M 
trials of Monte Carlo simulation, the conditional probability density can be written as 

o ^ 117 o ^ \ g{Ei,iii,ti,xi\Eo,no,to,xo) 
p{Ei,fti,ti,xi\Eo,no,to,xo) = ofjp r> . ^ > (Al) 

where g denotes the number of events that have an initial condition of Xq = (Eq, fl^, to, Xq) and a 
final condition of Xi = {Ei,Qi,ti,xi), and /3 denotes a distribution function of initial conditions 
of the simulation, which is normalized as 



f3{Xo)dXo = M. (A2) 

Then, we can write a distribution function of emitting photons in the system as 

/(Xi) = I p{Xi\Xo)a{Xo)dXo = I g{Xi\Xo)^^^dXo, (A3) 

where a{Xo) is a distribution function of the initial photons that are emitted by the source, or simply 
the source function. The simulation output is a list of the escaping photons with information on 
the initial and final conditions. For the list, f{Xi) is calculated by 

M 

f{Xi) = Y,m{X[,Xo), (A4) 

1=1 

{ "(-^o) (X' = x\ 
w,{X[,X^) = { P{Xo) ^ 1 (A5) 

I {X[^ X,) 

Here, Wi is a weight of each event and the summation is performed over all Monte Carlo events. 
When the geometry of the system is static or has a symmetry of time translation, it is possible to 
fix the initial time to in the simulation to zero; then Eq. (IA3P is reduced to 



fiEi,n,,h,x,)= [ giEi,ni,ti - to, xi\Eo, »o, 0, xo) "^"^ dEod^odtodxa. (A6) 



To simulate observations, we must consider time differences of light propagation to an observer 
from different positions. A distant observer directed toward an unit vector r observes the emission 
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at tobs = ti — xi ■ r/c, where c is speed of light and the origin of tobs is determined so that an 
emission at ti = and xi = is observed at tobs = through direct propagation (without any 
interactions). Events are selected to obtain observational spectra or images if the direction agrees 
with r within a certain small tolerance. The weight of the selected event is calculated by Eq. (jASP 
in which ti = tobs + Xi ■ r /c is substituted. If the system has spatial symmetry such as spherical or 
axial symmetry, we can extend the range of the direction selection to reduce statistical uncertainties 
unless the symmetry breaks. 



A. 2. Framework Design 

Our simulation framework, MONACO, provides methods for the Monte Carlo simulation and 
for the convolution to generate observed spectra or images ( ^A.ip . An outline of the calculation 
to obtain observational data is shown in Fig. [T3j The first part uses Monte Carlo simulations to 
calculate the probability of obtaining certain final conditions, given certain initial conditions. This 
part requires information about the initial conditions of the photon and the geometry including 
physical conditions of matter as input. The output here is a list entry for an event that has 
information about the initial state and the final state (emission). If a photon does not experience 
an interaction with matter, the final state is identical to the initial state. 

Using the event list generated by the simulation part as an input, the second part, referred 
to "convolution" in Fig. [T3l generates simulated observational data such as spectra or images by 
performing the convolution with the source function described in ^A.ll The output here is the final 
results, which can be compared with real observations. The list input can be a combination of 
different event lists to reduce statistical uncertainties or to extend the ranges of initial conditions. 
To perform the convolution by Eq. (IA5p . this part requires as input the initial conditions of photons 
both for the simulation and for the real source (source function). It also requires information on 
the observer; direction, time, distance and miscellaneous parameters. 

The first part (simulation part) performs the photon-tracking simulation including treatment 
of X-ray reprocessii ig. This part is an ex t ended and improved version of original Monte Carlo 
codes developed bv IWatanabe et al. (120031 . l2006l ) . We use the Geant4 simulation toolkit library 



(|AgostineUi et al.ll200,i 



Allison et al. 



20061 ) only for treatment of complicated geometry while phys- 
ical processes are our original implementations. The simulation framework consists of a manager 
and four user-defined modules. The manager, which is provided by the Geant4 library, handles 
communications between the modules and treats particle tracking based on Monte Carlo simula- 
tions and geometry navigation. The user-defined modules are for geometry, physical processes, 
primary generation, and data acquisition. The geometry module describes the geometry of the 
astrophysical object to be simulated as well as physical conditions of the matter. The module for 
physical processes organizes X-ray interactions with matter. The primary generator describes how 
to create initial photons to be tracked. For example, a point-like source is specified by a position 
vector of the source and an energy spectrum such as a power law with a photon index and an energy 
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Fig. 13. — Outline of our simulation framework 



range; in each Monte Carlo trial, the direction and energy of the primary photon are sampled from 
an isotropic distribution and the power-law distribution, respectively. The data acquisition module 
forms an event list of the simulation data which contains information on the initial and final states 
of the photon and saves it to a file. 



B. Scattering by Bound Electrons in Neutral Matter 



Sunvaev Churazovl (|l996l ) give a detailed review of scattering by bound electrons in neutral 



matter such as molecular clouds. Here we summarize physical implementation of this scattering 
process. This process can be divided into three channels according to the final state of the target 
electron: 

7i + Xj (ground state — Rayleigh scattering) 

7i + X* (excited state — Raman scattering) (Bl) 

7i + Xj_i + e~ (free state — Compton scattering) 

Assume that a photon is scattered by an electron bound by an atom. Through this scattering, 
the quantum state of the electron changes from an initial state |i) to a final state |/). The doubly 



7 + Xi 
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differential cross sections for atomic hydrogen are given by, 



hui + Ei = hi'2 + Ef, A/ii/ = hui — hi'2, 
q = ki-k2, x = q/h- 



(B2) 

(B3) 
(B4) 



Here, Ei and Ef are the initial and the final energies of the electron; hvi and ei are the energy 
and unit vector of polarization of the initial photon, respectively; hv2 and 62 are those of the final 
photon. The momentum transfer through this scattering is denoted by q; ki and ^2 are the initial 
and the final momenta of the photon. 

For atomic hydrogen, we have analytic solutions of the electron wave function. Thus, the 
differential cross sections of the three channels can be calculated by analytic formulae. For Rayleigh 
scattering, the differential cross section is 
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(B5) 



(B6) 



a is r-^/fi, and tb = h/nieca is the Bohr radius. For Raman scattering to the final electron state of 
n (principal quantum number), the differential cross section is given by 
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Substituting n = 1 to this expression of course yields Eq. (|B5p . 



Compton scattering by a bound electron makes a free electron, and the energy of the scattered 
photon is not determined uniquely by the scattering angle. The doubly differential cross section is 
given by 

d'^a / da\ 



dQ,dhv2 V / Th V ^^2 



1^1 l^]i?/.2 
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where 
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X [(gV + l-pV)2 + 4pV]-', 
//2m = - \Eb\ + /\hv. 
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Here, p and are the momentum of the ejected electron and the binding energy of the electron, 
respectively. The energy difference of the photon should be larger than the binding energy = 13.6 
eV for Compton scattering. 

Figure [H] shows the differential cross sections for the three channels of a 6.4 keV photon, 
corresponding to the iron Ka line. For small scattering angles, Rayleigh scattering is important, 
and Raman and Compton scattering are suppressed since the electron could not gain enough energy 
to be ionized or excited from the incident photon. On the other hand, Compton scattering becomes 
dominant at the large scattering angles and agrees with Compton scattering by a free rest electron 
described by Klein-Nishina formula. 




30 60 90 120 150 180 
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Fig. 14. — Differential cross sections for the three scattering channels of a 6.4 keV photon. The 
thin lines show contributions of Rayleigh scattering (red) , Raman scattering (green) and Compton 
scattering (blue). The thick black line is the sum of the total contributions. The values of the 
differential cross sections are normalized by the classical electron radius squared . 



In cold, dense matter, hydrogen exists in forms of molecules H2, not as atoms H. Since the elec- 
tron state of the hydrogen molecule differs from that of atomic hydrogen, it is important to consider 
the molecula r effects on the scatte ring processes. This problem under astrophysical applications is 
discussed by ISunvaev et al.l (jl999l ) . The most significant difference between a single-electron atom 
and a molecule is the enhancement of Rayleigh scattering due to electron coherence. Inelastic scat- 
tering (Raman and Compton) does not differ very much from that of atomic hydrogen. Moreover, 
effects of rotational and vibrational states of the molecule are negligible for t he X-ray energy band. 
Thus, we adopt a good approximation for molecular hydrogen proposed by ISunvaev et al.l (|l999l ) 
that scattering by the molecule is identical to that by atomic hydrogen except that the cross section 
of Rayleigh scattering is enhanced by a factor of two per electron. 



In many astrophysical situations, the second most abundant element, helium, also makes an 
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important contribution to scattering. Like molecular hydrogen, atomic helium has two electrons, 
resulting in enhancement of the coherent scattering. In addition, helium has larger binding energy 
-Eb = 24.6 eV and larger momentum of bound electrons than those of hydrogen. These affect the 
energy profile of scattered photons via Compton scattering. Since helium does not have analytic 
ex pressions of atomi c prop erties, we use numerical calculations of differential cross sections provided 
bv lVainshtein et al.l (11993). 



Figure [15] shows calculated spectra of scattered photons of a monochromatic 6.4 keV line via 
three different channels. Rayleigh scattering does not change the photon energy. There is a gap 
with a width of -^10. 2 eV that corresponds to the first excitation energy of hydrogen. Compton 
scattering is permitted for the scattered energy lower than E'une — = 6.4 — 0.0136 = 6.3864 keV 
due to hydrogen's binding effects. A discontinuous edge at 6.375 keV is due to the same effect by 
helium. This edge can be used for measurement of helium abundance for future large-area missions 
(jVainshtein et al.lll998l ) . The profile of the Compton scattering is smeared by the initial momentum 
of a target electron bound to an atom, displaying a long low-energy tail. If the electron were free 
and had zero momentum, the profile would have a sharp edge ai E = 6.24keV, corresponding to 
back scattering {9 = 180°). 
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Fig. 15. — Spectra of scattered photons of a monochromatic 6.4 keV line via Rayleigh scattering 
(red), Raman scattering (green), and Compton scattering (blue) in arbitrary units. 
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